Segregation in granular binary mixtures: Thermal diffusion 
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A solution of the inelastic Boltzmann equation that applies for strong dissipation and takes into 
account non-equipartition of energy is used to derive an explicit expression for the thermal diffusion 
factor. This parameter provides a criterion for segregation that involves all the parameters of the 
granular binary mixture (composition, masses, sizes, and coefficients of restitution). The present 
work is consistent with recent experimental results and extends previous results obtained in the 
intruder limit case. 
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Segregation and mixing of dissimilar grains is perhaps 
one of the most interesting problems in agitated granular 
mixtures. In some processes it is a desired and useful ef- 
fect to separate particles of different types, while in other 
situations it is undesired and can be difficult to control. 
Several mechanisms for segregation corresponding to dif- 
ferent scenarios have been proposed , but the problem 
is not completely understood yet. Among the different 
competing mechanisms, thermal diffusion becomes one 
of the most relevant mechanisms at large shaking ampli- 
tude where the sample of macroscopic grains resembles 
a granular gas. In this regime, binary collisions prevail 
and kinetic theory can be quite useful to analyze the 
physical mechanisms involved in segregation processes. 
Very recently, Schrdter et al. Q have carried out ex- 
periments in agitated mixtures constituted by particles 
of the same density. To the best of my knowledge, this 
is one of the few experiments in which thermal diffusion 
has been isolated from the remaining segregation mech- 
anisms. When convection is practically supressed, they 
report results to assess the influence of composition and 
dissipation on thermal diffusion. Given that this effect 
can be predicted by kinetic theory, in this Letter I study 
the particle segregation problem driven by the presence of 
a thermal gradient in order to explain some of the trends 
observed in the above experiments Q at large shaking 
amplitudes. This is the main motivation of this Letter. 

Under the above conditions, the so-called thermal dif- 
fusion factor Ay characterizes the amount of segregation 
parallel to the temperature gradient. Thermal diffusion 
is caused by the relative motion of the components of a 
mixture due to the presence of a temperature gradient. 
Due to this motion, concentration gradients subsequently 
appear in the mixture producing diffusion that tends to 
oppose those gradients. A steady state is finally achieved 
in which the separation effect arising from thermal dif- 
fusion is compensated by the diffusion effect. Here, I 
determine Ay from a solution [3| of the inelastic Boltz- 
mann equation that applies for strong dissipation and 
takes into account non-equipartition of energy. This lat- 
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ter feature, not present in the elastic case at equilibrium, 
has significant consequences on the transport coefficients 
of the granular mixture 0. 

The model system considered is a low density bi- 
nary mixture composed by smooth inelastic hard spheres 
(d = 3) or disks (d = 2) of masses rm and diameters <7i. 
Without loss of generality, we assume that ai > 02 . The 
inelasticity of collisions among all pairs is characterized 
by three independent constant coefficients of restitution 
c^iij 0^227 and CV12 = «2i- The mixture is in presence 
of the gravitational field g = —ge z , where g is a pos- 
itive constant and e z is the unit vector in the positive 
direction of the z axis. In experiments 0, the energy 
is usually supplied by vibrating horizontal walls so that 
the system reaches a steady state whose properties are 
presumed to be, far from the walls, insensitive to the 
details of the driving forces. Due to the technical dif- 
ficulties involved in incorporating oscillating boundary 
conditions, here particles are assumed to be heated by 
the action of a stochastic-driving force which mimics a 
thermal bath. Although the relation between this driven 
idealized method with the use of locally driven wall forces 
is not completely understood, it must be remarked that 
in the case of boundary conditions corresponding to a 
sawtooth vibration of one wall the condition to deter- 
mine the temperature ratio coincides with the one de- 
rived from the stochastic force |^| . This suggests that the 
results obtained from this driving method are a plausi- 
ble first approximation for qualitative comparisons with 
experimental results. As will be shown below, the re- 
sults derived here for the temperature ratio confirm this 
expectation. 

The thermal diffusion factor Ay (i ^ j) is defined at 
the steady state in which the mass fluxes j; vanish. Un- 
der these conditions, the factor Ay is given through the 
relation 

-Ay V In T = — V Xl , Ay + A# = 0, ( I ) 

where Xi — m/n is the mole fraction of species i (m is 
the number density of species i and n = n\ + 712). The 
physical meaning of Ay can be described by considering a 
granular binary mixture held between plates at different 
temperatures T and T' under gravity. For the sake of con- 
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FIG. 1: Temperature ratio T\jTi versus the size ratio a\ja-z 
for ay = a — 0.78 in the case of mixtures constituted by 
particles of the same mass density and equal total volumes of 
large and small particles. The lines are the kinetic theory re- 
sults in (a) the stochastic-driving case and (b) the free cooling 
case while the points refer to MD simulations Q. 
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FIG. 2: Phase diagram for BNE/RBNE for the case studied 
in Ref. [2j, namely, mixtures constituted by particles of the 
same mass density and equal total volumes of large and small 
particles. The data points represent the simulation results 
for a — 0.78 when convection is suppressed. Points below 
(above) the curve correspond to RBNE (BNE). 



creteness, we will assume that gravity and thermal gradi- 
ent point in parallel directions, i.e., the bottom is hotter 
than the top. In the steady state (ji = J2 = 0), Eq. 
describes how the thermal field is related to the composi- 
tion of the mixture. Assuming that A12 is constant over 
the relevant ranges of temperature and composition, in- 
tegration of Eq. (JU yields ln(xix' 2 /x 2 x[) = A i2 ln(T' /T). 
Here, Xi refers to the mole fraction of species i at the plate 
with temperature T and x\ refers to the mole fraction of 
species i at the plate with temperature T". Consequently, 
if T' > T and A l2 > 0, then x[ < x\ while if A l2 < 0, 
then Xi > x\. In summary, when A12 > 0, the larger par- 
ticles accumulate at the top of the sample (cold plate), 
while if A12 < 0, the larger particles accumulate at the 
bottom of the sample (hot plate). The former situation 
is referred to as the Brazil-nut effect (BNE) while the 
latter is called the reverse Brazil-nut effect (RBNE). 

The RBNE was observed first by Hong et al. |f| in 



MD simulations of vertically vibrated systems. They pro- 
posed a very simple segregation criterion that was later 
confirmed by Jenkins and Yoon [a by using kinetic the- 
ory. More recently, Breu et al. Q have experimentally 
investigated conditions under which the large particles 
sink to the bottom and claim that their experiments con- 
firm the theory of Hong et al. provided a number of 
conditions are chosen carefully. In addition to the verti- 
cally vibrated systems, some works have also focussed in 
the last few years on horizontally driven systems show- 
ing some similarities to the BNE and its reverse form . 
However, it is important to note that the criterion given 
in Ref. Jjj is based on some key assumptions: particles 
are assumed to be elastic, homogenous temperature and 
energy equipartition. These conditions preclude a com- 
parison of the present theory with the above simulations. 
Previous theoretical attempts to assess the influence of 
non-equipartition on segregation have been recently pub- 
lished. Thus, Trujillo et al. [9j have derived an evolution 
equation for the relative velocity of the intruders start- 
ing from the kinetic theory proposed Jenkins and Yoon 
that applies for weak dissipation. They use constitu- 
tive relations for partial pressures that take into account 
the breakdown of energy equipartition between the two 
species. However, the influence of temperature gradients 
which exist in the vibro-fluidized regime is neglected in 
Ref. because they assume that the pressure and tem- 
perature are constant in the absence of the intruder. A 
more refined theory has been recently provided by Brey 
et al. |10| in the case of a single intruder in a vibrated 
granular mixture under gravity. The present work cov- 
ers some of the aspects not accounted for in the previous 
theories |(J 0, since it is based on a kinetic theory 
that goes beyond the quasi-elastic limit 0,0 and applies 
for arbitrary composition x\ (and so, it reduces to the 
results obtained in Ref. 01 when x\ — ► 0). This allows 
one to assess the influence of composition and dissipation 
on thermal diffusion in bi-disperse granular gases with- 
out any restriction on the parameter space of the system. 
This the main value added of this paper and can be rele- 
vant to make comparisons with experiments/simulations 
in the dilute regime. 

To determine the dependence of the coefficient A12 on 
the parameters of the system, we focus our attention on 
an inhomogeneous steady state with zero mass flux and 
gradients only along the vertical direction (2 axis). Since 
the flow velocity vanishes, the momentum balance equa- 
tion yields d z p = —pg, where p = nT is the pressure, 
p = p± + P2, and pi = mini is the mass density of species 
i. To first order in the gradients, the constitutive equa- 
tion for the mass flux ji jZ is 

m 1 m 2 n o p , 

Ji.z = Dd z xx D p d z p - —D d z T, (2) 

P P T 

where D is the mutual diffusion coefficient, D p is the 
pressure diffusion coefficient, and D' is the thermal dif- 
fusion coefficient. The condition j± tZ — yields d z x\ = 
(p 3 /m 1 m 2 np)(D p /D)g - (p 2 /mim 2 p){D' / D)d z T '. Sub- 



3 



2.5 
2.0 

.{f 1.5 
1.0 
0.5 
0.0 







BNE 










(a) 


RBNE 





3 



FIG. 3: Phase diagram for BNE/RBNE in three dimensions 
for ctij = 0.7 and three values of composition: (a) xi = 0, (b) 
x\ — 0.3, and (c) x\ = 0.7. Points below (above) each curve 
correspond to RBNE (BNE). 



stitution of this relation into Eq. leads to 
np 2 D' - D p g* 



A 12 = 



P1P2 



D 



(3) 



where g* = pg/nd z T < is the reduced gravity acceler- 
ation. 

Since the mutual diffusion coefficient D is positive 
the sign of A12 is determined by the sign of the quantity 
D' — D p g*. Explicit expressions for the coefficients D' 
and Dp have been recently obtained from the Chapman- 
Enskog method in the first Sonine approximation 
Given that the driving-stochastic term does not play a 
neutral role in the transport, it must be remarked that 
the expressions for the transport coefficients obtained in 
the driven case slightly differs from the ones derived in 
the free cooling case. In particular, while in the former 
situation the coefficient D' = 0, the thermal diffusion 
coefficient D' = —(C,/2v)D p in the free case Q. Since 
the collision frequency v and the total cooling rate C are 
positive, then sgn(D') = sgn(— Dp). Consequently, the 
sign of A12 is the same as that of the pressure diffu- 
sion coefficient D p in the driving-stochastic case while 
sgn(Ai 2 ) = sgn(— D p {((/2v) + g*}) in the free case. 

Because of sgn(Ai2) = sgn(D p ) in the driven case, we 
will focussed henceforth on the explicit dependence of D p 
on the parameters of the system. Its expression is 



D p = 



PiP 



X 2 



p z v x 2 + £17 V/i 



1-1 



(4) 



where /i = m\jm 2 is the mass ratio, 7 = Ti/T 2 is 
the temperature ratio, and v is an effective collision frc- 
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The condition A12 = (or equivalently, D p = in 
the driven case) provides the criterion for the transition 
from BNE to RBNE. Equation (@J shows that the sign 
of Dp is determined by the value of the control param- 
eter 9 = "f/p.. This parameter gives the mean-square 
velocity of the large particles relative to that of the small 
particles. Thus, when 9 > 1 (9 < I), the thermal dif- 
fusion factor is positive (negative) which leads to BNE 



(RBNE). The criterion for the transition condition from 
BNE to RBNE is 7 = fi. In the case of equal granular 
temperatures (energy equipartition) , 9 — > \i~ x and so, 
segregation is predicted for particles that differ in mass, 
no matter what their diameters may be. It must be re- 
marked that, due to the lack of energy equipartition, the 
condition 9 — f is rather complicated since it involves all 
the parameter space of the system. In particular, even 
when the species differ only by their respective coeffi- 
cients of restitution they also segregate when subject to 
a temperature gradient. It is a novel pure effect of in- 
elasticity on segregation. The criterion for the transition 
BNE^=>RBNE is the same as the one fou nd p reviously 
in Ref . Q when otij is close to I and in Ref. |l0j | in the in- 
truder limit case (xi — > 0). However, as said before, the 
results obtained here are more general since they cover 
all the range of the parameter space of the system. 

To get the explicit dependence of the control parameter 
9 on otij , one still needs to compute the temperature ratio 
7. When the system is driven by means of a stochastic 
thermostat, the temperature ratio is determined from the 
condition |4| 



7C1 = pC,2 (white-noise thermostat), 



(5) 



where the explicit form of Q can be found in Ref. [3| . The 
condition JSJ) differs from the one derived in the undriven 
(free cooling) case where 7 is obtained by requiring 
that the partial cooling rates Q must be equal, i.e., 



Ci = C2 (free cooling). 



(6) 



This latter condition was used in Ref. [T(j to determine 
the difference of temperatures between an impurity and 
the surrounding gas in an open vibrated granular system. 

As said above, experiments and molecular dynamics 
(MD) simulations have been carried out very recently 
[2j to analyze segregation in agitated binary granular 
mixtures composed by spheres of the same material and 
therefore, the same mass density (i.e., /1 = (ci/o^) 3 ) 
and mechanical properties (oiy = a) . Figure ^ shows 
the comparison of the temperature ratio J1/T2 between 
MD simulations |2| and kinetic theory results based on 
the conditions JSJ) (driven case) and (undriven case). 
The experimental value of the coefficient of normal resti- 
tution is a = 0.78 and equal volumes of large and small 
particles are taken, i.e., x 2 = (01 /o^) 3 ^- While a good 
agreement between kinetic theory and MD simulations is 
found when the gas is assumed to be driven by a stochas- 
tic thermostat, significant discrepancies appear in the un- 
driven case, especially as the size ratio ai/a 2 increases. 
These results contrast with the ones obtained in Ref. 01 
in the tracer limit (xi — > 0) where the predictions of 
T\ /T 2 from kinetic theory based on the condition (f = Q 2 
compare quite well with MD simulations. However, it 
must noticed that for the cases studied in Ref. 0] the 
conditions (jSJ) and JSJ yield quite similar results for the 
dependence of the temperature ratio on the parameters 
of the system. 
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Consider next size segregation driven by thermal difu- 
sion. To make some contact with the recent experimental 
results of Schroter et al. [2| , let us consider again three- 
dimensional (d — 3) mixtures constituted by particles 
with the same mass density and equal total volumes of 
large and small particles. Figure [3 shows the phase dia- 
gram for this kind of systems. The results show that, 
for a given value of the coefficient of restitution, the 
RBNE is dominant at small diameter ratios. However, 
since non-equipartition grows with increasing diameter 
ratio, the system shows a crossover to BNE at sufficiently 
large diameter ratios. This behavior agrees qualitatively 
well with the results reported in Ref. 0] at large shaking 
amplitudes where thermal diffusion becomes the relevant 
segregation mechanism. At a quantitative level, we ob- 
serve that the results are also consistent with the sim- 
ulation results reported in 2J when periodic boundary 
conditions are used to suppress convection since they do 
not observe a change back to BNE for diameter ratios 
up to 3 (see red squares in Fig. 11 of Q). Although the 
parameter range explored in MD simulations is smaller 
than the one analyzed here, one is tempted to extrapolate 
the simulation data presented in Ref. Q to roughly pre- 
dict the transition value of the diameter ratio at a — 0.78 
(which is the value of the coefficient of restitution consid- 
ered in the simulations). Thus, if one extrapolates from 
the simulation data at the diameter ratios of 2 and 3, 
one sees that the transition from RBNE to BNE might 
be around o\jo2 = 10, which would quantitatively agree 
with the results reported in Fig. [21 Figure El also shows 
that the BNE is completely destroyed in the quasielastic 
limit (a ~ 1). 

Let us now investigate the influence of composition on 
segregation. Figure |31 shows a typical phase diagram in 
the three-dimensional case for = a = 0.7 and three 
different values of the mole fraction x\. The lines sepa- 
rate the regimes between BNE and RBNE. We observe 
that the composition of the mixture has significant ef- 
fects in reducing the BNE as the concentration of larger 
particles increases. In addition, for a given value of com- 
position, the transition from BNE to RBNE may occur 
following two paths: (i) along the constant mass ratio 



mxjmi with increasing size ratio ci/o^, and (ii) along 
the constant size ratio with increasing mass ratio Tfi^j vn\. 
Another aim in this paper is to assess the influence of dis- 
sipation on the phase diagrams. Our results show that in 
general the role played by inelasticity is quite important 
since the regime of RBNE increases significantly with dis- 
sipation. 

In summary, thermal diffusion (which is the relevant 
segregation mechanism in agitated granular mixtures at 
large shaking amplitudes) has been analyzed in the con- 
text of the inelastic Boltzmann equation by using a ki- 
netic theory that is not restricted to small dissipation and 
accounts for the non-equipartition of energy. The results 
reported here have been mainly motivated by recent ex- 
perimental results 0] focussed on the analysis of three dif- 
ferent segregation mechanisms: void filling, convection, 
and thermal diffusion. Concerning the latter effect, our 
model is able to explain some of the experimental and/or 
MD segregation results observed within the range of pa- 
rameter space explored. In addition, the theoretical pre- 
dictions for the temperature ratio of the mixture obtained 
from the driven condition (J5J agree very well with com- 
puter simulations. A more quantitative comparison with 
MD simulations in the low-density regime is needed to 
show the relevance of the present theory. As said before, 
comparison with MD simulations in the tracer limit case 
[x\ — » 0) 01 has shown the reliability of the inelastic 
Boltzmann equation to describe segregation. Given that 
the results derived here extends the description made in 
Ref. [10| to arbitrary values of composition, one expects 
that such good agreement is also maintained for finite 
values of X\. In this context, it is hoped that the present 
results stimulates the performance of such computer sim- 
ulations in the dilute regime. 
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